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ABSTRACT 

Sodium Laser Guide Stars (LGSs) are elongated sources due to the thickness and the 
finite distance of the sodium layer. The fluctuations of the sodium layer altitude and 
atom density profile induce errors on centroid measurements of elongated spots, and 
generate spurious optical aberrations in closed-loop adaptive optics (AO) systems. 
According to an analytical model and experimental results obtained with the Univer- 
sity of Victoria LGS bench demonstrator, one of the main origins of these aberrations, 
referred to as LGS aberrations, is not the Centre~of-Gravity (CoG) algorithm itself, 
but the thresholding applied on the pixels of the image prior to computing the spot 
centroids. A new thresholding method, termed "radial thresholding" , is presented here, 
cancelling out most of the LGS aberrations without altering the centroid measurement 
accuracy. 

Key words: instrumentation: adaptive optics - methods: analytical, laboratory. 



1 INTRODUCTION 

Sodium laser guide star (LGS) adaptive optics (AO) sys- 
tems allow in theory a full sky coverage; however, there are 
several limitations. The artificial star is elongated due to 
the sodium layer thickness (about 10 km) and finite dis- 
tance (90 km). Consequently, if the laser is launched from 
the secondary mirror holder, the spots of a Shack-Hartmann 
wavefront sensor (SH-WFS) are radially elongated from the 
centre to the edge of the pupil (Fig. [TJ. The spot elonga- 
tion is proportional to the telescope diameter and can reach 
several arcseconds for extremely large telescopes (ELTs). 

Moreover, the sodium layer is not static, b ut fluctuates 
with a time scale of about 1 minute or less (jPavis et al.l 
I2OO6I ). Fluctuations of the sodium layer altitude and atom 
density profile induce errors on centroid measurements of 
elongated spots and generate spurious aberrations on the 
wavefront in closed-loop AO systems. These aberrations, 
referred to as LGS aberrations, can reach several hundred 
nanometres peak-to-valley {ptv) for ELTs with the classical 
centre-of-gravity (CoG) c entroiding algorithm (jClare et al.l 
120071 : iLardiire et al.ll2008h . Some authors proposed new so- 
phisticated centroiding algorithms to mitigate LGS aberra- 
tions, such as the Matched Filtering (|Gilles and Ellerbroeh 
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|2008l: IConan et al.l 120091) or the correlation (|Povneerl I2OO3I : 
Thomas et al.l l2006. 2008). 

However, thanks to the LGS-bench demonstrator 
built at the University of Victoria (UVic) for ELT 
projects (|Lardiere et ahluOOSl : IConan ilaLlllooi), we found 
out that one of the main sources of the LGS aberrations 
was not the CoG algorithm itself, but simply the thresh- 
old applied on the pixels of the SH-WFS images before the 
centroid computation. 

Section [2] reviews the known possible origins of LGS 
aberrations. A model of the aberrations generated by thresh- 
olding is presented in Sec. |31 as well as a simple new thresh- 
olding method, termed "radial thresholding", which miti- 
gates most of LGS aberrations. The experimental results 
obtained with the UVic bench with the radial thresholding 
show that the CoG algorithm is still well suited for LGS 
wavefront sensing on ELTs (Sec. |4]). 



2 ORIGINS OF LGS ABERRATIONS 

If the LGS spots are radially elongated from the pupil cen- 
tre, as shown on Fig.[TJ the sodium layer fiuctuations induce 
centro-symmetric aberrations, such as focus (Z4) and spher- 
ical aberrations (Zn, Z22, etc.), and also square symmetric 
aberrations, such as tetrafoils (^14, ^26, etc.). 

The focus is due to a variation of the sodium layer aiti- 
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Figure 1. 29x29 elongated spot SH-WFS frame obtained with 
the UVic LGS-bench at SNR=18. The spot sampling is 2x8 pixel 
on the edge pupil, corresponding to 1" x4" on the sky, as expected 
for the TMT SH-WFS with a 1" seeing {aRON=0A3 ADU, no 
threshold). 



tude. This error is not an artefact and must be compensated 
by refocussing the LGS on the WFS with zoom optics and by 
updating the offsets of t he LGS WFS with a natural-guide- 
star (NGS) focus sensor ijHerriot et al . 2006). Non -common 
path errors of the LGS optical train, including the zoom op- 
tics, can vary with the sodium layer distance, i.e. the zenithal 
angle, and induce variable aberrations on the science path 
too. We assume that these systematic aberrations can be 
calibrated and virtually negated. 

Consequently, aberrations beyond focus are mainly 
artef acts of the wavefro nt sensing. According to a model 
from IClare et all (|2007l ) and to the first experimental re- 
sults obtained with the UVic LGS bench jLardiere et al.l 
|2008| ). the spherical aberrations arise due to a truncation of 
asymmetric LGS spots by a circular field-stop, while square 
symmetric aberrations are likely due to: 

• a spot truncation by a square field-stop or by pixel 
boundaries, horizontally and vertically elongated spots be- 
ing more truncated than diagonally elongated spots, 

• a spot overlap for square-grid lenslet arrays, horizon- 
tally and vertically elongated spots being more prone to 
overlapping, 

• quad-cell or sampling effects on centroid measurements. 

Both kinds of LGS aberrations have been reproduced 
and characterized in laboratory on the UVic bench with a 
time-series of 88 real sodium profiles (Fig. Beyond the 
focus, the most significative LGS aberrations detected are 
the spherical aberration Zn up to lOOnm ptv (30nm rms), 
and the tetrafoil Z14, with 40nm ptv (lOnm rms). Moreover, 
a correlation between the spherical aberratio n and the pro- 
file a symmetry was empirically established jLardiere et al.l 
|2008| ). Square-symmetric aberrations, such as Z14, should 
be mitigated by using a polar-coord inate CCD array 
l|Beletic et al.ll2005l : Ixhomas et aLll2008l ). 



However, we discovered later that Zn mode disappears 
if no threshold was applied on the pixels of the LGS WFS 
images before the computation of centroids. The threshold- 
ing discards the two extremities of each elongated spot, and 
consequently truncates radially each spot, as an optical cir- 
cular field-stop would do. The spot truncation caused by 
the field-stop is negligible compared to the truncation in- 
duced by the pixel thresholding if the field of view (FOV) 
of the outermost lenslets is wide enough to make an image 
of a 20 km-thick sodium profile. 

With such a large FOV and a polar-coordinate CCD 
array, as expecte d for the Thirty-Mete r-Telescope (TMT) 
LGS-AO facility (jEUerbroek et al.! '2008l . the thresholding is 
likely the main source of the LGS aberrations and deserves 
a specific study. 



3 MODELLING THE ABERRATIONS 
INDUCED BY THRESHOLDING 

Basically, a thresholding must be applied on the image pixels 
prior to computing the spot centroids in order to minimize 
the contribution of the detector read-out noise oron, or of 
the sky background. The thresholding is generally uniform 
over the pupil and is implemented as follow: 



h{x,y) = 



Thres 




I{x,y) ^ Thres 
I{x,y) < Thres 



(1) 



with / and It the raw and the thresholded images respec- 
tively. Thres is the intensity level of the threshold expressed 
in detector counts, i.e. in Analog-to-Digital Units (ADU). 
Generally, the threshold is defined from the readout noise (at 
3 a RON for instance). The remainder of this section demon- 
strates how a uniform thresholding can induce spherical 
aberrations when the spot profile is asymmetric. 

3.1 Geometrical model of the spot elongation 

If the sodium laser is launched on the system optical axis, 
the spots of a SH-WFS are radially elongated from the pupil 
centre, as shown on Fig. [21 Let p be the distance of a lenslet 
to the pupil centre, normaliz ed to 1 on the pupil edge. From 
Eq. 1 of lLardiere et al.l l|2008h . the geometrical radial size, on 
the detector plane, of the corresponding spot is proportional 
to p, such as E{p) ~ E p, with 



E ■ 



f RoNa 



/)2 
"0 



cos{z) . 



(2) 



In this equation, ho and ajva are the mean altitude 
and the thickness of the sodium layer respectively, z is the 
zenithal angle, R is the radius of the telescope pupil, and 
/ is the resultant focal length of the whole optical system. 
E is the geometrical radial extent, on the detector plane, of 
a spot located on the edge of the pupil (i.e. p = 1). E{p) 
is simply termed spot elongation hereafter, and can be ex- 
pressed in fim or in pixels. 

If X and y are the Cartesian coordinates of a spot in 
the detector plane, then the geometrical sizes of that spot, 
projected on the x and y axis are respectively defined as 
follow: 



Ex{x,y) = xE 
Eyix,y) = yE 



(3) 
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Figure 2. Geometrical model of the spot elongation if the laser 
is launched from the centre of the pupil. 



Let lo and a be the maximum intensity and the size of 
the central spot {i.e. p = 0) respectively. The extent a, ex- 
pressed in the same unit than E, is defined by the blurring 
due to the lenslet diffraction lobe, or to the seeing. Consid- 
ering this blurring effect, the actual length and width of a 
spot located at the radius p are respectively E{p) + a and 
a. Let Imax{p) be the maximum intensity of that spot. As 
the flux is conserved for all spots of the pupil, we have the 
following equation: 



^lo ~ Cr {E{p) + a) Imax(p) ■ 



Hence, 



(4) 



(5) 



3.2 Centroid errors versus threshold 

Let Po{r) be the projection of the sodium layer vertical pro- 
file on the SH-WFS detector plane, for a lenslet located 
on the edge of the pupil (i.e. p = 1). This projection is 
not linear and takes in account the perspective effect which 
stretches the lower part of the sodium layer and compress 
the upper part. Thus, the one-dimensional function Pair) 
is the radial intensity profile of the geometrical image of 
the most elongated spot. The variable r is the radius po- 
sition defined inside the sub-image plane of the considered 
spot only. If a threshold of value Thres is applied on the 
spot image, the thresholded profile, obtained accordingly to 
Eq.[TJ is noted Pt{r), with t = Thres / Imax{p) the threshold 
value normalized to the local spot maximum intensity. For 
any spot profile and any threshold t £ [0 ; 1] , we can define 
the dimension-less function g{t) such as 

1 // rPt(r)dr J r Po{r) dr 
^ ^ ^ V / Pt{r)dr / Po{r)dr 



(6) 



According to our previous assumption, this function 
provides, for any spot of the pupil, the radial centroid error 
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Figure 3. Radial profile of the theoretical geometrical image of 
the most elongated spot for a sodium profile featuring a strong 
asymmetry (profile #77 of the sequence displayed in Fig.|6j. This 
one-dimensional curve is the projection of the vertical sodium 
profile on the SH— WFS detector plane, sampled on a finite num- 
ber of pixels, with no diffraction, no atmospheric turbulence, no 
photon noise and no read-out noise. 



induced by the threshold t, expressed as a fraction of the 
elongation of the considered spot. In absolute units (pixel 
or pm), the radial centroid error scales as the spot elonga- 
tion, and is g{t)E{p). We note that g{t) — for symmetric 
profiles. Figure [3] plots a geometrical spot profile Po{r) with 
its function g{t). This profile is among the most asymmetric 
sodium profiles found in the 88-profile sequence of Fig. |6l 
The projection of the radial error g{t) E{p) on the x and 
y-axis, provides 5x{x,y) and 5y{x,y), the errors made on 
the wavefront slopes in x and i/-directions respectively for a 
spot located in coordinates (x;y): 



5x{x,y) = Ex{x,y)g{t) = x E g{t) 
Sy{x,y) = Ey{x,y)g{t) = y E g{t) 



(7) 



Hence, the wavefront error Z(x, y) caused by threshold- 
ing is deduced from the slopes, such as. 



(8) 



Using a polar coordinate system (p, 6), such as a; = 
p cos 9 and y = p sin d, we can use the following conversion 
formulae. 



P dp ~ ^ dx ~^ y dy 



I sz _ 
to obtain. 



y dx 



+ X 



dZ 

dy 



az 

dp 

az 
as 



pEg{t), 



= 0. 



(9) 



(10) 



Thus, the wavefront distortion generated by threshold- 
ing is always a surface of revolution whatever g(t), i.e. what- 
ever the sodium profile. In terms of Zernike polynomi- 
als (|Nolllll976D . this result proves that the thresholding can 
only induce a combination of focus {Z4) and spherical aber- 
rations {Zii, Z22, Zsr, etc.). 

3.3 Wavefront aberrations induced by uniform 
thresholding 

If the threshold value Thres is constant for the whole pupil, 
then t — Thres /Imax (p) varies with p. From Eqs. [5] and [TOl 
we obtain in that case: 
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dZ ^ (Thres / E 

Moreover, we assume th; 
wise approximated, as a polj 



15- 



5W 



git 



with the sum starting at i = 
tion (Eq. [SJ . Using the binon 
integration of Eq. 1111 



n 



Thres \ ' 



lo 



J A^ j + 2j\{i^jy. 



(3) 



(13) 



We can check that Z{p) = if the spots are not elon- 
gated {E = 0), or if the threshold is null, or again if the 
spot profile is symmetric {gi = 0,Vi). If n is the order of 
the polynomial g{t), the radial order of the wavefront aber- 
ration is n -I- 2. A uniform threshold apphed on a triangular 
or trapezoidal (i.e. n = 1) asymmetric spot profile induces 
already aberrations of the third radial order: 



^, , E Thres 



2E 



(14) 



The same threshold applied on a more complex asym- 
metric sodium profile will induce higher order aberrations 
on the wavefront. The projection of on the Zernike poly- 
nomials base is a pure focus: 

P^^^Z., (15) 

while p'^ is a combination of focus and all spherical aberra- 
tions (Zii, Z22, Zz7, etc.): 
3 



0.3 {Zi + 0.14 Zu - 0.015 Z22 + 0.003 Z37 + . 



(16) 



The focus term is corrected by the NGS focus sensor al- 
ready required to offset the altitude variation of the sodium 
layer. Beyond the focus, the main aberration induced by 
a uniform thresholding is Zu mode, higher order spheri- 
cal aberrations being negligible. The Zn error is fluctuating 
with the time, due to the variations of the sodium profile 
asymmetry, and is an issue for LGS AO systems. 



3.4 Wavefront aberrations induced by radial 
thresholding 

If the threshold value is now, not uniform over the whole 
pupil, but defined for each sub-image proportionally to the 
maximum intensity of the local spot, then the normalized 
threshold t = Thres / 1 max [p) becomes constant, as well as 
g{t). From Eq. [10] we deduce that 



Z{p) = ^g{t)p 



(17) 



Consequently, such a thresholding method induces only 
a focus error and no spherical aberrations at all, whatever 
the sodium profile structure. As the threshold value scales 
as the inverse of the radius, this thresholding is refereed to 
as "radial thresholding". 
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Figure 4. Horizontal radial cut of the 29x29 elongated spot 
SH-WFS frame of Fig. [T] obtained with the UVic LGS-bench 
at SNR=18. The red solid line illustrates a uniform threshold set 
at 3 times the read-out noise (a-/jojv=0.43 ADU on the bench). 
The radial cut is a binning of the 3 central pixel lines of the frame. 
According to the spot elongation, the theoretical intensity ratio 
between the central spots and the edge spots is 4. The actual 
ratio is around 7 here due to the gaussian profile of the coUi- 
mated beam used on the bench, and also to the field curvature of 
the imaging lens. These efl'ects exaggerate slightly the wavefront 
errors induced by thresholding. 
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Figure 5. Truncation effect induced by thresholding on the 4 
outer spots at SNR=18, for 3 difl'ercnt threshold levels. The 0- 
ADU level is the detector offset. Pixel values might be negative 
due to read— out noise. 



4 LABORATORY RESULTS OF THE RADIAL 
THRESHOLDING 

4.1 Experimental setup 

The radial thresholding is easy to implement on real SH- 
WFSs, and has been validated in laboratory with the UVic 
LGS demonstrator bench and compared to the results ob- 
tained without threshold or with an uniform threshold. Two 
signal-to-noise ratios (SNR) have been considered, 140 and 
18, the latter being expected for real LGSs. At SNR=140, 
the read-out noise is negligible, and thresholding could be 
not necessary. However, this favourable case is useful to 
properly distinguish the thresholding efi^ect from the noise, 
and to validate our model. 

The readout noise of the SH-WFS camera is 
0.43 ADU rms. As the offset frame of the detector is sub- 
tracted to the raw images, some pixels may have negative 
values due to the read-out noise. Consequently, a 0-ADU 
threshold discards the negative values, and does not mean 
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Na profile sequence (data from TMT.AOS.TEC. 07.005 Nap) 



Uniform threshold, SNR=140 
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Figure 6. Sodium profile sequence used on tlie UVic LGS— bench 
(real data from Purple Crow LIDAR, Univ. of Western Ontario, 
altitude X time resolution is 25 mx70 s) and LGS aberrations 
measured by a centre— of— gravity— based SH-WFS in open-loop 
with no turbulence, when no threshold, or a uniform threshold 
(UTCoG), or again a radial threshold (RTCoG) is applied on the 
pixels of LGS spot images. The focus {Z4) component due to an 
altitude change of the sodium layer is removed here to highlight 
the component due to the thresholding (SNR=140, UT=13 ADU, 
RT=4% of the local spot maximum, crjjQjY=0.43 ADU). 



that no threshold is applied on the image. We will have to 
distinguish both cases in our following experimental study 
at low SNR. 

The FOV of the circular field-stop of the bench corre- 
sponds to a sodium layer thickness of 20 km for the most 
elongated spot. As the mean thickness of the sodium layer is 
about 10 km, the spot truncation effect due to the field-stop 
should be negligible compared to the spot truncation caused 
by the thresholding. 

Figure[T]displays a 29x29 radially elongated spot image 
obtained with the UVic bench for SNR=18. On this image, 
the maximum spot elongation _E is 8 pixels, and the spot 
width o" is 2 pixels. Figure |4] plots an horizontal cut of this 
image {y = 0). A uniform threshold at 1.3 ADU, i.e. SaaoN 
for the bench, is shown as an example. This graph already 
points up that a uniform thresholding is not suited for LGS 
spots, since edge spots are fainter than the central spots. At 
low SNR, a uniform threshold might even reach the maxi- 
mum intensity of some spots located on the edge. Figure [5] 
illustrates this issue for SNR=18. 
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Figure 7. LGS-induced aberrations measured on the UVic LGS— 
bench at SNR=140 when a uniform threshold (a) or a radial 
threshold (b) is applied on LGS spot images (a-j;ojv=0.43 ADU). 



4.2 Performance at high flux 

Figure [6] plots the LGS aberrations for Z4,, Zn, Z22 and Z14 
modes versus the time, measured during the whole sodium 
profile sequence, with no threshold, with a 13-ADU uniform 
threshold and a radial threshold (4% of local maximum) at 
SNR=140, in open-loop without turbulence. Such a high 
uniform threshold, not suited for real AO systems, is chosen 
here to exaggerate the thresholding effects on the wavefront 
and validate the model. Figure [7] plots the corresponding 
temporal RMS fluctuation of all Zernike modes beyond focus 
up to Z2S, for SNR=140 and for different threshold values 
(both uniform and radial). 

As the main focus error source is the altitude variation 
of the sodium layer, the focus obtained without thresholding 
has been subtracted from those obtained with thresholding. 
This subtraction highlights the extra focus error induced by 
thresholding. This error reaches 350nm ptv with the uni- 
form threshold, and 200nm ptv with the radial threshold. 
Moreover, a correlation is visible between the focus and Zn 
modes induced by the uniform thresholding. This correla- 
tion confirms that Z4, and Zn modes arise together. This 
result is in good agreement with the model (Egs. [T3lto ll6|l . 
Although the ratio Z^/Zn depends on the coefficients gi, 
I.e. the sodium profile structure and asymmetry, we can at- 
tempt to quantify it for a linear fit of the functions g{t), 
I.e. for a sequence of trapezoidal or triangular sodium pro- 
files. In this simple case, the ratio Z4/Z11 expected by the 
model is about 7.3 with E/a = 4 (Eqs. [l4]to[l6}, while the 
observed ratio is about 4.9 for the ptv amplitudes (Fig. |6]). 
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Uniform threshold, SNR=18 
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Radial threshold, SNR=18 
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Figure 8. LGS-induced aberrations measured on tlie UVic 
LGS-bench at SNR=18 when a uniform threshold (a) or a 
radial threshold (b) is applied prior to computing centroids 
(cTRON=0-i3 ADU). 



The model and the experiments provide the same order of 
magnitude for the ratio Zi/Zu. The difference is Ukely due 
to the assumption made about the profile shape. A higher 
order fit of the functions g(t) would be required for each 
profile, as well as further measurements using a symmetric 
sodium profile as a reference. 

Figure [6] shows clearly that the radial thresholding miti- 
gates the spherical aberrations {Zu and Z22) induced by the 
sodium layer variations, as if no thresholding was applied. 
The variations for the Zu mode reaches 90nm ptv with uni- 
form thresholding, against lOnm ptv with radial threshold- 
ing or without thresholding at all. The remaining error is 
likely the component induced by the noise and by the circu- 
lar field-stop. Moreover, Fig.[7K shows that the RMS fluctu- 
ation of Zii mode scales as the level of the uniform thresh- 
old, which is in agreement with the model too (Eg. I13|) . 

As expected, the radial thresholding does not mitigate 
the square symmetric LGS aberration Z14. However, the Z14 
error is modified by the thresholding, and is almost null 
without thresholding (Fig. [6] and (T^)- Hence, the threshold- 
ing would likely induce other effects not considered in our 
model. The pixel boundaries or the spot overlap may have 
more or less impact depending on the threshold value. As- 
suming that polar coordinate CCD arrays should mitigate 
square symmetric modes, we will not investigate with more 
details the other possible sources of this mode in this paper. 

At SNR=140, the best performances are obtained with 
a radial threshold set at 4% of the local spot maximum. The 



Table 1. RMS wavefront errors (WFE), in nm rms, and achiev- 
able Strehl ratios (S), in % for A = l^m, obtained at SNR=140 
and 18, during the sodium profile sequence with the center- 
of-gravity algorithm without pixel thresholding, with uniform 
thresholding, and with radial thresholding. The threshold levels 
are optimized for SNR=18. The uniform threshold values are 13 
and ADU for SNR=140 and 18 respectively, while the radial 
threshold values are 4% and 12% for SNR=140 and 18 respec- 
tively. All Zernike modes beyond the focus (from Z5 to Z2s) are 
added up and spatially averaged over the pupil. The bench ac- 
curacy is the ultimate centroiding performance achieved by the 
bench on non-elongated spots. Readout noise is 0.43 ADU rms. 





SNR= 


140 


SNR= 


18 




WFE 


S 


WFE 


S 


No threshold 


4.34 


99.9 


60.87 


86.4 


Uniform threshold 


10.65 


99.5 


31.91 


96.1 


Radial threshold 


4.27 


99.9 


18.70 


98.6 


Bench accuracy 


2.73 


99.9 


17.84 


98.7 



RMS wavefront errors averaged over the pupil for all modes 
beyond the focus are summarized in Table [1] 



4.3 Performance at low flux 

Figure |8] plots the temporal RMS fluctuation of Zernike 
modes from Z5 to Z2S, for SNR=18 and for different thresh- 
old values, both uniform and radial. 

At low SNR, edge spots are very faint, a few ADU. Con- 
sequently a uniform threshold, even low, can easily exceed 
the maximum of some spots, and induce huge errors. The 
errors induced by a 3-ADU threshold exceed 100 nm, the 
upper limits of the graph (Fig.|H^). Hence, a small change of 
the threshold level or of the SNR, which is fluctuating with 
LGS, can decrease significantly the performances. A 0-ADU 
uniform threshold appears to be more robust in that case, 
but it is not optimal in term of noise. 

The performance obtained with the radial threshold is 
more robust at SNR=18, and does not depend much on the 
threshold value (Fig.|8]D). This is certainly due to the ability 
of the radial thresholding to track the temporal variations of 
the spot intensity: the radial thresholding is also a dynamic 
thresholding. The wavefront errors measured with a radial 
threshold are lower than those measured with a uniform 
threshold in the best case (Tab.[T}. At SNR=18, the optimal 
value for the radial threshold is 12% of the local maximum. 

Lastly, Tab.[T]shows that the performance achieved with 
the radial thresholding on elongated spots is very close to 
the ultimate performance achievable on the bench with a 
non-elongated source. This result basically means that the 
radial thresholding cancels out most of the errors induced 
by the LGS spot elongation, without altering the centroid 
accuracy in presence of noise. 



5 DISCUSSIONS 

5.1 Weighted Centre— of— Gravity algorithm 

Considering the issues caused by the thresholding on LGS 
spots, we could suggest to use the Weighted Centre-of- 
Gravity algorithm, which does not require a prior thresh- 
olding, instead of the conventional centre-of-gravity algo- 
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Z(p) = + 8/3 LGS aberrations with uniform thresholding 




Zernike modes 



Figure 9. Typical LGS aberrations arising wlien tlie laser is 
launched from the side of the pupil (x=0 and y=0) with a uniform 
thresholding (top) and for a radial thresholding (bottom) . These 
aberrations are compared, in term of Zernike modes, to those aris- 
ing when the laser is centrally projected. Aberration amplitudes 
are peak— to- valley in an arbitrary unit. 



rithm. Basically, this algorithm consists in weighting the 
pixels of the spot ima ge by a reference image of the 
spot (jNicoUe et al.|[20o3 ). If the spot image matches per- 
fectly the reference image, this technique is equivalent to 
compute the center-of-gravity of the square of the image 
intensity. However, the centre of gravity of the square of an 
elongated spot image differs from the actual center of gravity 
of that spot, if the spot profile is asymmetric. Consequently, 
the Weighted-Center-of-Gravity algorithm cannot be used 
on LGS AO systems as it is. 



5.2 Case study of non— radial spot elongation 

The theoretical model presented in Sec. [3] assumes that 
the laser is launched on the optical axis. This configura- 
tion is optimal for ELTs because it minimizes the spot 
elongation and the LGS aberrations. However, on a IO- 
meter class telescope, the laser can preferably be launched 
from t he side of the telescope pupil, like for the Keck tele- 
scope (jWizinowich et al.|[200^ . The LGS spots of a SH- 
WFS are not longer radially elongated from the pupil centre, 
but from an edge of the pupil. 

The model is still valid for this case, providing that 
the origin of the coordinate system {p,0) or {x,y) coincides 
with the position of laser launch telescope. Equation 1131 
expressing the wavefront Z{p,9), remains unchanged and 
the aberrations generated by a uniform thresholding are still 
centro-symmetric for a virtual giant pupil centred on the 
laser axis, with a diameter twice the distance between the 
laser axis and the optical axis. The only difi'erence is the 
domain of definition of the function Z{p,8), which must be 
shifted to match the telescope pupil. 

Figure |9] displays the portion of the wavefront Z(p, 6) 
seen by a circular pupil located at the right of the laser axis, 
and plots the peak-to-valley amplitudes of this wavefront 
projected on Zernike modes. Assuming a linear fit of the 
function g{t), we can consider Z{p) oc p^ + 8/3 for repro- 
ducing the LGS aberrations induced by an uniform thresh- 
olding (Eq. [14] with E/a = 4), and Z{p) oc p^ for reproduc- 



ing the LGS aberrations induced by a radial thresholding 
(Eq.[T7j. 

As the aberrations are not centro-symmetric, the 
Zernike decomposition of p^ features, in addition to the fo- 
cus and spherical aberrations, many other modes, such as tip 
and tilt {Z2, Z3), astigmatism (Z5, Zq), coma (Z7, Zs) and 
some trefoil (Zg, Ziq). Moreover, compared to the centrally 
projected laser configuration, the peak-to-valley amplitude 
of these aberrations are about 7 times greater. 

On the other hand, the decomposition of p^ features 
only tip and tilt modes in addition to the focus, meaning 
that the radial thresholding removes all the LGS aberrations 
beyond the focus. Tip-tilt and focus modes are not issues, 
since they will be corrected by the NGS tip-tilt and focus 
sensor already necessary to determine the atmospheric tip- 
tilt, unseen by a LGS WFS, and to offset the sodium layer 
distance variations, respectively. 

Finally, the radial thresholding is even more advised if 
the laser is launched from the side of the pupil, since the LGS 
aberrations induced by the uniform thresholding are much 
higher and are spread in more modes. The side-projected 
laser configuration has no been tested on the UVic bench, 
but could be implemented if the need arises for any LGS AO 
project. 



6 CONCLUSION 

The fluctuations of the sodium layer vertical profile in- 
duce centro-symmetric and square symmetric aberrations 
on LGS SH-WFSs. Centro-symmetric aberrations, such as 
focus and spherical aberrations, are due to a spot trunca- 
tion by either a circular fleld-stop or by a threshold applied 
on the image pixels. Square symmetric aberrations, such as 
Zi4, are due to a spot truncation by the pixel boundaries 
and a sampling effect. Some extra aberrations may arise due 
to uncalibrated non-common path errors of the LGS train 
varying with the sodium layer distance. 

Both kind of aberrations have been reproduced in the 
laboratory with UVic LGS bench. Square symmetric aber- 
rations should be negated with a polar coordinate CCD 
array. Considering such a detector and a field of view per 
sub-aperture wide enough to be able to image a 20-km 
thick sodium layer, the residual LGS aberrations are mainly 
caused by the threshold applied on the pixels of elongated 
spots prior to computing the centroids. 

This statement is confirmed by the theoretical model 
and experimental results presented in this paper. The aber- 
rations induced by the thresholding contain focus and spher- 
ical aberration modes if the laser is launched from the pupil 
centre. Their amplitudes scale as the spot elongation {i.e. the 
telescope diameter), the sodium profile asymmetry and the 
threshold value. 

These aberrations disappear if the threshold value is 
not static and uniform over the pupil, but defined dynam- 
ically and independently for each lenslet, proportionally to 
the maximum intensity of the local spot. The residual LGS 
aberrations are only a focus error if the laser is centrally 
projected, or a combination of tip, tilt and focus if the laser 
is launched beside the pupil. These residual modes are not 
issues since they will be corrected by the natural guide star 
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tip-tilt and focus sensor already required in any LGS AO Wizinowich P., Le Mignant D., Bouchez A., et al., 2006, 

system. PASP, 118, 297 

This new thresholding method, termed radial thresh- 
olding, is very simple and add no extra computational time 
compared to the uniform thresholding. It can be imple- 
mented on all current LGS SH-WFS in operation. More- 
over, the radial thresholding makes the centre-of-gravity 
algorithm still well-suited for Extremely Large Telescopes, 
where the LGS spot elongation is greater. Unlike the 
matched filter or the correlation centroiding algorithms, the 
centre-of-gravity is a simple, well-proven and fast algorithm 
which requires no special calibration processes or reference 
images, and has no repercussions on the design of the LGS 
AO system. 
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